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(57) Abstract 

A method fo imaging within a sample space an object containing 
quadrupolar nuclei comprises irradiating the object to excite nuclear 
quadrupole resonance, applying to the object a magnetic field gradient 
having a profile (B 0 (x)) such that the square of the profile ((B 0 (x))2) va- 
ries linearly with distance (x), detecting resonance response signals 
from the nuclei, and deriving an image from the response signals (see 
Figures 1 to 7). An analogous apparatus is also disclosed. A method 
of and apparatus for nuclear quadrupole resonance testing an object, 
and a method of and apparatus for detecting the presence of a particu- 
lar substance containing a given species of quadrupolar nucleus, are 
also disclosed (see Figures 7 and 8). 
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METHODS AMI) APPARA TUS FOB HOB TFCTTHfi 

The present invention relates to a method of and apparatus 
for imaging within a sample space an object containing 
quadrupolar nuclei. The object may, for Instance, be a human or 
animal body or a part thereof; the image may be a simple 
spin-density image or may be a temperature, pressure or electric 
field distribution Image. For convenience the imaging referred 
to above will be called NQR imaging. The present invention also 
relates to a method of and apparatus for NQR testing an object, 
and to a method of and apparatus for detecting the presence of a 
particular substance containing a given species of quadrupolar 
nucleus. 

Methods for producing Images from the resonances of spin-y, 
nuclei which have a magnetic moment but no electric quadrupole 
moment have been extensively developed and described 1n various 
books (e.g. P. Mansfield and P.G. Morris. "NMR Imaging In 
Blomedlcine", 1982, Academic Press). Such methods may use linear 
magnetic field gradients superimposed on a strong, uniform 
20 magnetic field. The resonance frequencies of the spin-fc nuclei 
are linearly dependent on the total magnetic field; this is made 
to vary linearly with distance across the sample, and the 
distribution of the substance concerned is derived from the 
frequencies in the signals produced. Such techniques are 
25 commonly called NMR imaging. 

Methods for Imaging using the resonances of quadrupolar 
nuclei (which have 1*1) would be expected to have several 
advantages. Firstly, the quadrupolar resonances can be detected 
without using the strong, uniform magnetic field which 1s needed 
to make the magnetic resonances of spin-X nuclei conveniently 
detectable. Hence the relatively large, expensive and 
sample-size-limiting magnet structures which are necessary for 
NMR imaging will not be needed. Avoiding the need for this 
strong magnetic field also avoids the substantial complications 
35 and possible distortions caused by non-uniformity In the field 
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and magnetic interference produced in practical NMR equipment. 
Secondly, the quadrupolar resonances are more definitely or 
characteristically associated with specific chemical 
environments, and it is therefore easier to distinguish results 
5 which are due to a particular substance from effects which are 
due to other substances present in the sample. Methods In which 
quadrupolar resonances might be observed directly rather than 
through their interactions with magnetic resonances would be 
highly suitable for medical use because they do not require 

10 patients to be exposed to rapidly changing strong magnetic fields. 

Possible disadvantages or limitations on NQR imaging are that 
the nuclei may be less abundant, may have nuclear quadrupole 
resonance frequencies lower than the magnetic resonance 
frequencies of the protons commonly used in NMR imaging, and may 

15 have smaller gyromagnetic ratios. This can create a sensitivity 
problem which has two aspects. Firstly, it may be difficult to 
achieve a degree of excitation of the quadrupolar nuclei 
comparable with the excitation of the spin-% nuclei- which is 
commonly used in NMR equipment. The radiofrequency power used in 

20 medical applications at least must be limited to avoid undue 
heating and damage to living tissue. Secondly, the response 
signals will be weak and of low signal-to-no1se ratio, needing 
sophisticated data processing for their detection. 

To our knowledge methods for imaging using quadrupolar 

25 resonances have as yet had very little development. Rommel £± al 
(J. Magnetic Resonance 21. 630-636, 1991) reported theoretical 
reasoning with a conclusion that there would be great practical 
difficulties in any attempt to use such methods to derive any 
image from the interaction of a magnetic field gradient with 

30 quadrupolar nuclei, considering in particular half-integral spin 
nuclei of I£3/2. As an alternative they reported a method in 
which the strength of a radiofrequency excitation was varied 
across a sample in zero magnetic field. 

Matsui fit lL <J- Magnetic Resonance £&. 186-191, 1990) 

35 reported measurements on 35 C1 nuclei in sodium chlorate in a 
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linear magnetic field, in which an image was derived from the 
effects of the field on the spectral shape and width of the 
resonance. It appears that this may be a special case; the 35 C1 
nuclei have I«3/2 and their asymmetry parameter 1n sodium 
5 chlorate Is zero. Matsui i± &L state that "it is practically 
impossible" for the magnetic field to shift the resonance 
frequency. This may be incorrect, at least in relation to nuclei 
other than 35 C1 in substances other than sodium chlorate. 

According to one aspect of the present invention, in general 

10 terms there are provided methods and apparatus for producing an 
image representing the distribution of a specific substance or 
local temperatures within an article or in part of a human or an 
animal body, wherein the article or part of the body is 
irradiated with a pulse or pulses of radlofrequency energy at or 

15 close to a quadrupole resonance frequency of atomic nuclei of I£l 
in the said specific substance and is subjected to at least one 
magnetic field gradient having a profile which produces a 
variation (including possibly a splitting) of the said resonance 
frequency which is a linear function of distance in the direction 

20 of the said magnetic field gradient, response signals from the 
said nuclei are measured and an image is derived by data 
processing the result of many such tests. According to the 
invention, preferably either the magnetic field profile is a 
function B 0 (x> of distance x such that its square CB 0 (x)] 2 1s a 

25 linear function of x, within a limited sample volume of the 
article or body, or else the invention relates to a method of 
imaging using quadrupolar nuclei having integral spin quantum 
number (I o ], 2, 3, ...), preferably having unity spin quantum 
number. 14 N is a nucleus for which I - 1. 

30 According to an aspect of this invention, there is provided 

apparatus for imaging within a sample space an object containing 
quadrupolar nuclei, comprising means for irradiating the object 
to excite nuclear quadrupole resonance, means for applying to the 
object a magnetic field gradient having a profile B such that B n 

35 varies linearly with distance, n being greater than one, means 
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for detecting resonance response signals from the nuclei; and 
means for deriving an image from the response signals. An 
analogous method is defined in Claim 1. 

This invention arises from our discovery that, surprisingly, 
5 for a wide range of substances containing quadrupolar nuclei, a 
gradient profile having n greater than one (rather than a linear 
gradient) should be used in NQR imaging since it (rather than a 
linear gradient) produces a variation in resonance frequency 
which is a linear function of distance. As will be apparent from 
10 the following, for many (possibly the majority of) applications 
we have discovered that n should be 2 (i.e. a square law 
gradient), by which is meant that n should be exactly, or close 
to, two. 

Magnetic field gradients in different directions may be 

15 applied in sequence to derive sufficient results for the 
derivation of a two-dimensional or three-dimensional image, each 
gradient having a profile which makes the resonant frequency vary 
as a linear function of distance in the direction of the 
gradient. Where gradients in only two directions are used, the 

20 sample or article may be rotated about one of these directions to 
derive information sufficient for three-dimensional imaging. 

The sensitivity and excitation difficulties noted above may 
in practical terms restrict the species of nuclei which can be 
imaged, and limit the size of article or body parts to which the 

25 present technique may be applied. The present invention is 
particularly concerned with producing images of the distribution 
of various specific substances containing 14 N nuclei. Such 
nuclei are present in collagen, urea, peptides, proteins, amino 
acids and many other substances which are important in normal or 

30 abnormal human anatomy and/or metabolism, and images representing 
their distribution in body parts may be useful in the diagnosis 
or assessment of disease and in monitoring medical treatments. 
It also appears feasible to use the methods and apparatus of this 
invention to produce Images of substances containing 23 Na, 39 K, 

35 35 C1 or 127 I nuclei. For other nuclei such as 2 H, 7 L1, 9 Be, 11 B, 
17 0, 25 Mg, 27 A1, 51 V, 55 Mn, 59 Co, 75 As, 79 Br, 81 Br, 197 Au and 
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209 B1, the sensitivity and excitation problems may be more revere 
and with presently available techniques it is possible that any 
imaging of these nuclei would be limited to comparatively small 
and/or inanimate samples. 
5 The quadrupole resonance frequencies are also shifted by 

temperature variations, and by comparing tests with different 
magnetic field profiles the technique may be extended to produce 
images representing the variations of temperature in a sample or 
body part. This can provide a non-invasive method for measuring 
10 or monitoring localised temperatures within a sample or body 
part, and it may be useful in monitoring and/or controlling 
treatments which involve localised heating. Imaging of pressure 
or electric field variations may be effected following the same 
principle. 

15 In the sub-molecular environment of compounds or crystals, 

the nature and disposition of the adjacent electrons and atomic 
nuclei produce electric field gradients which modify the energy 
levels and the resonance frequencies of the quadrupolar nuclei. 
The symmetry or asymmetry of these local electric field gradients 

20 is an important factor in any study of nuclear quadrupole 
resonance, and it is measured by an asymmetry parameter n which 
can have values from zero to one. For instance, for the ring 
nitrogen in the explosive RDX, n*0.62, whilst for the peptide 
nitrogen in amino acids and polypeptides n*0.4; the nitrogen in 

25 cocaine and heroin is thought to have r^O.l. The spin quantum 
number I is another controlling factor. 14 N and 2 H nuclei have 
1-1; nuclei such as 1A B and 75 As have I«3/2, and 17 0 and 27 A1 
1*5/2, etc. 

The case of nuclei having integral spin quantum number is 
30 considered first. 1=1 nuclei are taken as a representative 
example, but the principles behind the following analysis, and 
the conclusions drawn from the analysis, apply to all integral 
spin quantum number nuclei. In the case of nuclei having I«l 
there will be three resonant frequencies « x , u y and « z which will 
35 always have a relationship 

<o x =»y+« z . ( i ) 
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When first and second-order perturbation theory are applied to a 
study of these resonance frequencies in an externally applied 
magnetic field, it appears that there are substantially two cases 
which can be represented by different mathematical approximations. 
5 For such nuclei in substances having n>0.1, over a 

considerable range of magnetic field strengths the resonance 
frequencies of individual nuclei may be expressed by equations, 
derived from a paper by S. Pissanetzky (Zeeman effect of 14 N 
nuclear quadrupole resonance with polycrystalline samples, J. 

10 Chem. Phys., Vol. 59, No. 8. p. 4197, 1973), of the form 

« x = (3+n>K + (G 2 B 2 /K).F X (2) 
«v = (3-n>K + (G 2 B 2 /K).F y (3) 
4 - 2 n K + <G 2 B 2 /K).F Z (4) 
where G « y/2v* y = the gyromagnetic ratio, and K * |e 2 q0/4h|, 

15 e 2 qQ/h being the quadrupole coupling constant. F x , Fy and F 2 are 
functions involving n and the two polar angles (6,4>) specifying 
the relative orientation of the local electric field gradient 
with respect to the direction of the applied magnetic field B. 
Since in general there will be many molecules randomly oriented 

20 with respect to the direction of the applied magnetic field, 
there will be many nuclei with a range of values of F x , F y and F z 
and consequently a range of resonant frequencies from each 
location with any non-zero value of B; the resonance is 
broadened. However, the maximum value of the dominant term in 

25 the expressions for F x , F y and F z is proportional to 1/n, and 
hence this broadening is comparatively limited for substances 
which have a relatively high asymmetry parameter Ti- 
lt is possible with much computation to calculate how the 
resonances of a population of randomly-oriented nuclei will 

30 combine to determine the line shape and frequency shift of their 
observable resonance effect. 

The results of our researches are that, in general, the 
observable resonances will show a frequency shift proportional to 
B 2 with a gradual broadening of the line shape, until a critical 

35 field strength B 2 is reached; above this field strength the line 
shape becomes distorted and the frequency shift gradually departs 
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from a strict B 2 dependence. The angle a between the directions 
of the magnetic field B and the radiofrequency excitation field 
B x has a significant effect on the onset and nature of these 
changes, and it affects different lines differently. 

We have conducted experiments using pure substances in a 
uniform magnetic field of adjustable strength to study how the 
field strength alters the line shape and shift. From such 
experiments using samples of various explosives containing 14 N 
nuclei such as RDX (for the ring 14 N of which n«0.62>, we have 
found that the resonant frequency « x does shift in proportion to 
B 2 over a substantial range of magnetic field strengths. Me have 
also found that above a certain field strength the resonances we 
have studied develop a doublet structure when B is perpendicular 
to B x . Fine structure is also indicated in the results published 
by Matsui fit al. though they did not comment on it. In the case 
of RDX, both components of the doublet shift in proportion to B 2 , 
with different rates; we would expect the same to be true for 
other substances. 

Hence for NQR imaging of nuclei having 1-1 in substances 
having n>0.1, it is preferable to use a magnetic field profile 
B(x) such that 

EB(x)] 2 » u - vx (5) 
where u and v are numerical constants. This profile will clearly 
make the resonant frequency vary as a linear function of 
distance x. The same principle applies to nuclei with 1-2, 3,... 
etc., provided that the appropriate frequency is selected. 

We have discovered that such a profile may be produced by the 
combined effects of a central coil centred around a sample space 
and one or more offset coils coaxial with the central coil but 
offset to one side of the sample space, when the total 
ampere-turns product of the one or more offset coils is 
substantially greater than the ampere-turns -product of the 
central coil, and the fields of the offset coils and the central 
coil act in the same direction. 

We have also discovered an alternative arrangement which uses 
two coils of different diameters, both offset to one side of the 
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sample space, when the ampere-turns product of the larger coil Is 
substantially greater than the ampere-turns product of the 
smaller coil, and the field of the smaller coil opposes the field 
of the larger coil. 
5 Still considering the case of integral spin quantum number 

nuclei, substances having n zero or close to zero (say n<0.1) are 
now considered. Taking again the example of I«l , a different 
approximation leads to the following equations (see the 
Pissanetzky paper) 

io w x = 3K + ( n 2 K 2 + G 2 B 2 cos ?B>* (6) 

co y = 3K - ( n 2 K 2 + G 2 B 2 cos 2e)* (7) 
« z « 2 ( n 2 K 2 + G 2 B 2 cos 2 6)^ (8) 
for the resonance frequencies of individual nuclei. Assuming 
that the nuclei are randomly oriented, we have deduced that in 

15 this case the observable resonance peak will shift in proportion 
to the field strength B provided that n is small enough or, more 
particularly, provided that n 2 K 2 «G 2 B 2 cos 2 6 (that is, provided 
that (GBcos6/(nK)) 2 »l). 

However, for many substances of practical interest for which 

20 n<0.1, (GBcos0/(nK)) 2 can be actually considerably less than 
one. For instance, in a 10 Gauss field and assuming n«0.05, this 
ratio equals 0.002 for Heroin and 0.04 for Cocaine. When, as in 
these cases, this ratio is much less than one, then we have 
discovered, by use of the Binomial expansion, that, as with 

25 n>0.1, the resonance frequency will shift according to B 2 . The 
principles adduced above also apply to nuclei with 1-2, 3,... etc. 

Hence, for a majority of the substances of practical interest 
with nuclei having integral spin quantum number, a magnetic field 
profile of the form given in Equation 5 is appropriate. 

30 However, there may be some substances for which 

(GBcos0/(r|K)) 2 is near unity and for which therefore the above 
analysis does not hold. Sometimes, dependent on the substance, a 
profile substantially of the form given in Equation 5 would still 
be appropriate, with any slight deviations from a precise B 2 

35 dependence being accommodated by appropriate minor adjustments to 
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the field gradient profile to give a linear frequency shift. In 
other circumstances the deviations from a B 2 dependence may be so 
significant that a completely different gradient profile may need 
to be employed. Such a profile would normally be somewhere 
5 between square law and linear (that is, l<n<2, where n Is the 
power of the profile referred to previously), although there may 
be some cases where n 1s actually greater than two. 

The case of nuclei having half integral spin quantum number 
is now considered. I - 3/2 nuclei are taken as a representative . 
10 example, but the principles behind the following analysis, and 
the conclusions drawn from the analysis, apply to all half 
integral spin systems. Nuclei with 1-3/2 in a single crystal 
structure in a weak magnetic field B have four quadrupole 
resonance frequencies given by equations of the form 

15 "n = u o * EB - F n (e . n> (9) 

where E s a function of y and n and the four resonances have 
different functions F n . This leads for single crystals to 
frequency shifts linearly dependent on the magnetic field 
strength B. 

20 Hence a linear magnetic field profile will be appropriate for 
imaging nuclei with 1=1, 2, 3 ... etc. when (GBcose/( n K)) 2 »l or 
nuclei with 1-3/2, 5/2, etc. in a single crystal structure. This 
can be provided by coils in the anti-Helmholtz or Maxwell 
configuration - two colls of radius R, spaced R apart and 

25 connected to produce opposing fields. A system providing linear 
gradient fields in three orthogonal directions will have three 
pairs of Identical coils symmetrically arranged about the sample 
space. In such a system complementary magnetic field profiles 
can be provided by simply reversing the direction of the 

30 currents. Such systems may be suitable, for instance, for the 
imaging of parts of human arms or legs and possibly heads; they 
should be simpler and less expensive than present NMR Imaging 
instruments. 

For substances with I - 3/2, 5/2, etc. other than single 
35 crystal structures, the dependence of resonance frequency on 
magnetic field becomes complex and not readily amenable to 
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theoretical analysis. In these circumstances, we propose- that 
the relationship between resonance frequency and magnetic field 
is predetermined for a specific substance of interest and stored 
as a "look-up" table. Then a magnetic field profile of the form 

5 given in Equation 5 (or otherwise a linear field gradient) may be 
employed to produce spatial encoding. The image can be derived 
by appropriate use of the look-up table. Alternatively, a field 
profile such as would produce a variation of resonance frequency 
which is substantially a linear function of distance 1n the 

10 direction of the magnetic field gradient could be utilised. This 
profile would usually have a form somewhere between linear and 
square law (that is, l<n<2). 

By employing the magnetic field profiles described above 
known methods of NMR imaging can be adapted to NQR imaging, the 

15 major simplification being the elimination of the strong 
background magnetic field which is necessary for NMR experiments. 

In NMR imaging it is usual to use radiofrequency pulses of 
sufficient power to cause a 90° rotation of the magnetisation of 
the nuclei. The pulse duration t w is given by the equation 

20 * 

t w « (io> 

2 y B, 

where B x is the peak strength of the radio frequency excitation 
25 field. There will be other constraints on the pulse duration and 
so in effect this equation often sets the power required. In NQR 
experiments the corresponding condition for polycrystalline 
samples containing spin-1 nuclei is 

2ir 

3yB x 

The gyromagnetic ratio y will usually be much smaller for 
quadrupole nuclei than for protons, so that the r.f. field 
35 strength and power required to satisfy equation (11) will be much 
greater than that required to achieve the corresponding 
conditions in NMR imaging of protons. Hence it may be 
practically preferable or necessary to use pulses of considerably 



30 
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lower power and rotation angle, and to make up for this by 
accumulating results from more pulses, and probably applying 
pulses more frequently, so that the Interval between pulses (t) 
1s less than the sp1n-latt1ce relaxation time (T^. This may be 
5 particularly necessary to avoid tissue damage in medical 
applications. This will also make the accumulated signal 
strength less sensitive to variations In the rotation angle over 
the sample space. Shaped or phase/amplitude modulated pulses can 
also assist in reducing the r.f. power. 

10 The projection-reconstruction methods of Imaging developed 
from the work of Lauterbur (Nature 242, 190-191, 16/3/73) may be 
preferred for medical applications as they require only small 
step-wise changes in the magnetic field gradients. 

Variations of temperature may cause frequency shifts 

15 comparable with the frequency shifts caused by the magnetic 
field. This will cause distortion in images of samples which 
have an inhomogeneous temperature distribution, but this can be 
corrected. This effect may also be used to advantage as a 
non-invasive means for Imaging or monitoring the temperature 

20 distribution within the sample. Over a small temperature range, 
the variation in frequency is often linear with respect to 
temperature . 

Fourier transformation of a set of results taken in zero 
magnetic field with a sufficiently short r.f. pulse will 
25 therefore give a spectrum which, when compared with one or more 
spectra of a homogenous uniform temperature sample of the 
relevant substance, will Indicate the range of temperatures 
present. 

Increasing temperature generally shifts the resonances to a 
30 lower frequency, so that in any image produced by excitation of 
the frequency oj x without temperature correction, Indications of 
the nuclear density in any hot spot will appear displaced towards 
the region of lower magnetic field. Repeating the experiment 
with the field profile reversed (or complementary) will cause the 
35 hot-spot signal to appear displaced In the opposite direction, 
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and comparison with the previous image will indicate the position 
of the hot spot and its temperature. 

Variations of pressure or electric field will also cause 
frequency shifts which may be utilized for imaging purposes. 
5 Where the magnetic field splits any resonance into a 

raultiplet, or where there are two or more resonant frequencies 
close together, it is advantageous * to use excitation signals 
having a frequency distribution broad enough to excite all 
resonances or all parts of the raultiplet. The response signals 

10 will then show a pattern of slow beats at the difference 
frequencies. In certain cases, e.g. when the line is a simple 
doublet, it may be desirable to use the frequency of these slow 
beats as the characteristic to be measured, rather than the 
frequency of either resonance, because it will usually be more 

15 sensitive to variations in field strength and it may be more 
readily distinguishable. 

The disclosure of the preceding paragraph gives rise to a 
further aspect of the invention which is now discussed. This 
aspect concerns the method of and apparatus for NQR testing and 

20 the method of and apparatus for detecting the presence of a 
particular substance containing a given species of quadrupolar 
nucleus referred to previously. It has particular relevance to 
the detection of 14 N nuclei in explosives (such as RDX, HMX, TNT 
or PETN) or narcotics (such as Cocaine or Heroin). The invention 

25 might, for example, be embodied In a testing instrument installed 
at an airport. However, it could, in appropriate circumstances, 
find use in other areas, such as imaging. 

It is known (see T. Hirschfeld and S.M. Klainer, "short range 
remote NQR measurements", J. Molecular Structure, 58 (1980) 

30 63-77) to conduct NQR tests to detect the presence of a 
particular substance containing a given species of quadrupolar 
nucleus (such as 1A N nuclei) by exciting, using a repeating 
radiofrequency pulse sequence, nuclear quadrupole resonance at a 
single resonance frequency of the nuclear species and by 

35 detecting the response signals at that resonance frequency. 
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However, such testing techniques can be inaccurate and can 
lead either to the detection of too many "false positives" or to 
the detection of too few actual positives. 

According to the further aspect of the present invention, 
5 there is provided a method of nuclear quadrupole resonance 
testing an object, comprising repetitively irradiating the object 
to excite during each repetition nuclear quadrupole resonance at 
a plurality of resonance frequencies, and detecting the response 
signals at each of the resonance frequencies. 

10 By repetitively exciting and then detecting NQR resonance at 
a plurality of resonance frequencies the accuracy of the testing 
can be increased without there necessarily being a concomitant 
proportional increase in test time. The repetition time (t) of 
the irradiation 1s generally limited by the spin-lattice 

15 relaxation time (T x ) of the substance of interest; t may 
typically need to be at least 3T X to ensure full signal recovery 
between repetitions. A number of resonances can be excited and 
their response signals detected within this repetition- time. If 
these response signals are summed, signal-to-noise ratio will 

20 generally Increase according to the square root of the number of 
resonance frequencies excited. Thus the Invention can afford 
greater accuracy than has hitherto been achievable without any 
significant increase 1n test time. 

Suitably, the Irradiation may either take the form of 

25 radiofrequency pulses centred at a single frequency and having a 
frequency distribution sufficiently broad to excite all of the 
relevant resonances (e.g. shaped or phase and/or amplitude 
modulated pulses) or take the form of radiofrequency pulses at or 
close to each of the resonance frequencies of Interest. 

30 To ensure a rapid test, the repetitions preferably occur at 

least once every 5 T x (where T t is the spin-lattice relaxation 
time of the NQR nuclei of Interest) and usually more frequently 
than that (e.g. i - 3T X ). This can ensure that the repetitions 
are no longer than 1s necessary to achieve full signal recovery 

35 between repetitions. Where T 2 varies according to which 
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resonance is being excited, the terra "spin-lattice relaxation 
time" suitably connotes the longer or longest value of T x . 

One possible variant of the invention is that the testing is 
not repetitive, but that tests at the plurality of resonance 
5 frequencies overlap. 

Hence, preferably, irradiation at one of the resonance 
frequencies begins before detection at another of the frequencies 
ends. By overlapping the excitation and detection of the various 
resonances in this way a somewhat quicker test can be afforded 
10 without loss of accuracy. The excitation and detection of the 
various resonances may suitably be substantially simultaneous. 

Preferably, the resonance frequencies are well -resolved, so 
that there is little or no energy transfer between the various 
resonances and the maximum signal-to-noise ratio can be 
15 attained. However, even if this condition is not fulfilled, 
signal-to-noise ratio improvements can still be achieved to the 
extent that saturation of all the overlapping resonances does not 
occur. 

Preferably, the resonance frequencies are components of a 

20 multiplet. A multiplet usually connotes two or more resonance 
frequencies which are relatively close together (usually within a 
few percent of each other) and are yet we 11 -re solved. Multiplets 
may arise, for instance, due to crystal lographic non-equivalence 
of otherwise structurally identical NQR nuclei (as is the case 

25 with, for example, RDX, TNT and HMX); they may also arise (e.g. 
Heroin and Cocaine) due to the nearly axial ly symmetric electric 
field gradient in which the NQR nuclei are located. 

The particular advantage of applying the invention to 
substances in which the resonance frequencies are components of a 

30 multiplet is that the resonance frequencies are usually neither 
so far apart that more than one radiofrequency probe (coil) and 
pre-amplifier are required to excite and detect the various 
resonances nor so close that energy is transferred between one 
resonance and another. It will be understood that a probe's 

35 capability to irradiate at different frequencies depends on its 
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Q-factor. A typical probe might be able to irradiate over a 
range of 50kHz without retuning and 500kHz with retuning. 
Retuning can be effected In a fraction of a second and may 
therefore not adversely affect the duration of the tests. 
5 However, even when the resonance frequencies are not 

components of a multiplet it is still possible to make use of the 
present invention. For simplicity, the lowest of the resonance 
frequencies may suitably be, say, at least 501, preferably 751, 
of the highest. With these percentage limitations it may be 

10 possible to use the same rf probe and pre-amplifier (although 
suitable adjustments to both will possibly be necessary). As an 
example, the explosive PETN (which exhibits no multiplets) has 
resonance frequencies near 360 and 460kHz, the lower being within 
roughly 801 of the higher. 

15 According to a different facet of this further aspect of the 

invention, there is provided a method of detecting the presence 
of a particular substance containing a given species of 
quadrupolar nucleus, comprising Irradiating the object to excite 
nuclear quadrupole resonance at a plurality of resonance 

20 frequencies of the given nuclear species, detecting response 
signals at each of the resonance frequencies, and providing an 
alarm signal in dependence on whether nuclear quadrupole 
resonance is detected at each of the plurality of resonance 
frequencies. 

25 By this arrangement, the number of "false positives" (alarm 

signals erroneously indicating the presence of the substance) can 
be considerably reduced. If only one resonance is excited and 
detected, a false positive can readily occur since there 1s 
Inevitably some outside noise or Interference which can produce 

30 "spikes" or other features on the instrument output which may 
appear to be genuine NQR response signals at or near the relevant 
resonance frequency. However, if the alarm signal- is provided in 
dependence on whether a resonance of the appropriate intensity is 
detected at each of the plurality of resonance frequencies, the 

35 probability of a false positive occurring can be considerably 
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reduced, especially if testing takes place at a significant 
number of resonance frequencies (say, at least three or four). 

Preferably, the alarm signal is provided in further 
dependence on whether the relative intensities of the response 
5 signals at each of the plurality of resonance frequencies match 
those which would be expected for the particular substance. This 
measure can further reduce the number of false positives. It is 
envisaged that this measure may suitably be put Into practice by 
predetermining in separate tests the relevant characteristics of 

10 the particular substance and storing these for use 1n the field. 

In a preferred embodiment, the response signals from the 
different resonances are summed together, and if, on the basis of 
this summation, the presence of the substance 1s detected to a 
given (usually statistically predetermined) degree of 

15 probability, a second test Is performed In which further 
Irradiation repetitions are effected and the alarm signal 1s 
provided in dependence on whether resonance Is detected (again to 
a given, probably different, degree of probability) at each of 
the plurality of resonance frequencies. If the presence of the 

20 substance is not detected to the given degree of probability in 
the first test, the second test can be omitted. This two-stage 
procedure can thus afford a quick and yet accurate technique for 
detecting the presence of a particular substance. 

Since resonance frequency varies with temperature, pressure 

25 or magnetic field, it may necessary to take such variations into 
account by irradiating over a range sufficient to cover all 
possible variations of resonance frequency which might in 
practice be encountered. 

The further aspect of the invention extends to apparatus 

30 according to Claims 35 to 39. 

This further aspect of the Invention is now exemplified with 
reference to the NQR testing of specific substances. 

In the case of RDX (cyclo-trimethylenetrinitramine) the 
multiplicity produced by crystal lographic non-equivalence of 

35 otherwise structurally identical nitrogen nuclei can be 
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exploited. In the vapour phase, the molecule Is said to have C 3v 
symmetry, but this is lost in the solid state, and all three 14 N 
quadrupole resonance lines from the ring nitrogen nuclei are 
well-resolved triplets. The splittings are small (less than 150 
5 kHz for v x near 5,100 kHz), and, although the individual 
frequencies vary with temperature, the splittings themselves are 
almost independent of temperature between -30 and +30°C, varying 
from 143,52 kHz at the former to 145,49 kHz at the latter. 

As suggested above, the signal-to-noise ratio can be 

10 increased and the number of false positives considerably reduced 
by, for example, sequentially pulsing at these three frequencies, 
separately collecting the free inducation decays (f.i.d.'s) 
and/or echoes and processing each of them. Addition of the 
resulting Fourier transformed responses or their areas will then 

15 increase the signal-to-noise ratio by a factor of /n, where n is 
the number of lines in the multiplet. (It will be appreciated 
that, in this further aspect of the invention, "n" has a meaning 
different from that in the first aspect.) For RDX,- n-3. By 
interleaving pulses at each individual frequency, as is described 

20 later, the same time can be devoted to accumulating those three 
signal responses as would be needed for just one. Since the 
signal-to-noise ratio has been improved by a factor of /3, 
however, the statistical certainty of detection is greater, and 
the number of false positives can be reduced. 

.25 After each pulse or pulse sequence at a given frequency, the 

NQR -instrument will analyse the absorption peak for its 
frequency, height, width and area, taking into account possible 
errors assessed from the noise measured off -resonance. In this 
example, only if statistically equivalent responses are detected 

30 after each excitation frequency has been excited would the 
instrument provide an alarm signal. 

The above technique can be highly discriminatory against 
noise peaks and interference; only RDX will give equal, or 
nearly equal, response signals at the three selected frequencies 

35 at a given temperature. 

Heroin and Cocaine are now considered. We have discovered 
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that Heroin shows a close doublet (v x , v y ) near 3,900 kHz and 
Cocaine a doublet near 960 kHz. In these cases, the splitting 
arises from the nearly-axially symmetric electric field gradient 
in which the 14 N nucleus is located, and not from resolved 
5 crystal lographic non-equivalence. The present invention still 
applies, and the signal-to-noise ratio may be increased by /2 in 
these two cases, although the lines are not completely 
independent. 

TNT in both Its monoclinic and orthorhomblc forms gives six 
io v x lines between 850 and 900 kHz due to crystal lographic 
non-equivalence, so n-6. HMX gives two v x lines between 5000 and 
5300 kHz. 

We have discovered that PETN gives only one 14 N v x signal 
near 460 kHz, but v y is nearby at 360 kHz and it may be possible 
15 to excite both effectively in the same probe, in which case n«2. 

Embodiments of the invention will now be described, by way of 
example only, with reference to the accompanying drawings, of 
which:- 

Figures 1 and 2 are diagrammatic representations of 
20 alternative electromagnet coil configurations which may be used 
to produce a magnetic field profile B(x) such that [B(x>] 2 varies 
as a linear function of distance x; 

Figures 3 and 4 are respectively a front elevation and an end 
elevation of apparatus for NQR imaging; 
25 Figures 5 and 6 are graphical diagrams representing signals 

used in the operation of the apparatus; 

Figure 7 is a block circuit diagram of electrical equipment; 

and 

Figures 8 are timing diagrams for an embodiment of the 
30 invention. 

Figure 1 shows in axial section two coils 1 and 2 which have 
the same mean radius R and are spaced a distance R apart on a 
common axis 3. In use it will be arranged that coil 1 will have 
twice the ampere-turns (ni) of coil 2, the current directions 
35 being the same so that their fields are in the same direction. 
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If distance along the axis from the plane of coll 1 1s denoted by 
x, this arrangement produces a magnetic field profile B 01 (x> such 
that tB 01 <x)3 2 = bj-ajX, where b l and a t are constants, for 
values of x in the range 0.6R<x<1.3R; this desired profile Is 
5 provided within a cylindrical sample space indicated by the 
broken lines 4. Any nuclei with 1=1 which are placed within this 
space will have precession rates aird resonance frequencies « x 
which will vary as a linear function of the coordinate x, 
provided, where n <0.1, that the ratio (GBcose/( n K)) 2 Is much less 
10 than one. 

Figure 2 shows two concentric colls 6 and 7 which have mean 
radii of 2R and R respectively. In use it may be arranged that 
the ampere-turns (3n1> of the outer coil is three times the 
ampere-turns (ni) of the inner coil, and the field of the inner 

15 coil opposes the field of the outer coll. This arrangement 
produces a field profile B 02 (x) such that [B 02 (x>] 2 = b 2 - a 2 x 
where b 2 and a 2 are constants, for values of x in the range 
R<x<1.9R. The drsired profile is produced within a cylindrical 
space indicated by the broken lines 8 and again any nuclei in 

20 this volume with 1=1 will have precession rates and resonance 
frequencies u x which will vary as a linear function of the 
coordinate x, with the proviso mentioned in the previous 
paragraph. 

The arrangement of Figure 1 will clearly require much less 
25 electrical power to establish a given range of magnetic field 
strengths than the arrangement of Figure 2, but the fact that its 
sample space 4 lies mainly In the centre of a coil and does not 
extend very far from the coil may be Inconvenient for some 
applications. The Figure 2 arrangement has the advantage that 
30 Its useful sample space 8 is well separated from the coils and 
more accessible. This will make it more convenient for studies 
on human or animal patients. When the Figure 2 configuration is 
used it will be important to ensure that no part of the sample or 
body or any material which might have nuclear quadrupolar 
35 resonances 1n the frequency range being used should be present in 
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the space where x<R, between the colls and the sample space* 8, as 
this would produce confusing results in the same frequency range 
as the desired signal from the sample volume of the article or 
body under test. 

5 The exact spacing size and power relationship of these coils 

is not critical and both designs may be varied to make the sample 
space size and position suitable for particular applications or 
to minimise the power required; many similar arrangements may be 
used to achieve similar results. 

10 In both arrangements, slight adjustments to the number of 

ampere-turns and/or the relative disposition of the two coils can 
be made to accommodate slight deviations from precise B 2 
dependence which may occur in particular compounds. 

Figures 3 and 4 show a front elevation and an end elevation 

15 of apparatus for NQR Imaging using coils in the configurations of 
Figure 1 and Figure 2 to produce appropriate magnetic field 
gradients in three orthogonal directions. For the sake of 
clarity the end elevation Figure 4 is drawn to a larger scale 
than the front elevation Figure 3. 

20 A sample object 11 is placed on a cradle or platen 12 which 

runs on rollers 13 along a rail 14. The rail 14 is mounted on 
pedestals 15. An open-ended cylinder 16 of electrically 
insulating material surrounds a central portion of the rail so 
that the cradle 12 and sample objects 11 can pass through it. 

25 This serves to shield the sample objects from any malfunction 
such as arcing or fire which could occur from insulation failure 
or overheating in the electrical equipment, and it supports r.f. 
coils 17-20, with associated tuning components (not shown). 
These r.f. coils are saddle-shaped and mounted on the outside of 

30 the cylinder 16. They will be used to excite and to detect 
nuclear resonance signals. It may not be necessary to have two 
pairs of colls; one pair (17-19 or 18-20) may suffice or 
alternatively a tuned r.f. coil may be wound like a solenoid on 
the cylinder 16. The cradle 12, rail 14, cylinder 16 and other 

35 structural parts will be made of materials selected to avoid 
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distortion of the magnetic or r.f. field and to avoid confusing 
resonances. 

Two shielded electromagnet coils 21 and 22 are placed around 
the cylinder 16, to provide a magnetic field profile along the 
5 longitudinal . direction of the rail 14, substantially in 
accordance with the principle of Figure 1. The asymmetry of 
these coils should be noted; the smaller coll 22 is placed over 
the central part of cylinder 16 and the r.f. colls 17-20, and the 
larger coll 21 is displaced to one side by an amount equal to the . 
10 mean radius of the colls. 

Two concentric electromagnet colls 23 and 24 are supported 1n 
a horizontal plane, centred Just below the bottom of coll 22. 
The Inner coil 24 Is hidden by the outer coil 23. These colls 
will be used to provide a magnetic field profile 1n the vertical 
15 direction through the centre of the volume within the 
cylinder 16. Two similar concentric coils 25 and 26 are mounted 
on a support structure 27. in a vertical plane Immediately behind 
the coll 22. The inner coll 25 1s hidden by the outer coll 26 In 
Figure a and by the other components in Figure 3. Both these 
20 coil pairs are arranged and used substantially 1n accordance with 
the principles of Figure 2 but coil 23 is made a little smaller 
and coll 26 a little larger so that the rear of coll 23 can fit 
inside the lower part of coil 26. The colls 25 and 26 provide a 
magnetic field profile In a transverse horizontal direction. 
25 Positions within the sample space will be represented by 
orthogonal coordinates x, y and z in the directions indicated by 
axes Ox, Oy and Oz. 

Figure 5 shows how the magnet currents may be varied In part 
of a typical cycle of tests. The traces show the current 1 x 
passed through the coils 25 and 26, the current 1 y passed through 
coils 21 and 22, and the current 1 z passed through the colls 23 
and 24. The part of the cycle shown comprises six scanning 
periods in which the currents are held constant, separated by 
transition periods 1n which one current Is decreased while 
35 another Is increased. Each current varies from zero to a maximum 



30 
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value I. The annotations below the scanning periods indicate the 
directions of the field gradients produced, the notation xy 
indicating a gradient along the direction of the line x=y, z«0. 
The time scale of the whole cycle in medical applications will be 
5 determined by the maximum allowable rate of change of magnetic 
field to avoid any risk of adverse effects. The whole cycle may 
typically take from one to ten minutes and will include several 
hundred scan periods. 

During each scan period a pair of the r.f. coils will be 

10 excited with pulses at a repetition rate determined by the 
relaxation times of the quadrupole frequency to be detected and 
of relative phases determined by the nature of the pulse cycle. 
For example, one half of the excitation pulses in each scan 
period (preferably alternate pulses) may be phase-inverted. The 

15 resulting free induction decay and/or echo signals will be 
sampled and measured in sampling periods occurring at a set time 
after each excitation pulse and before the next excitation 
pulse. Preferably fast digital sampling is used to take many 
measurements in each sampling period. Measurements made at 

20 corresponding times are accumulated and averaged but the 
responses from the phase-inverted excitations are subtracted from 
the responses to the non-inverted excitations. Quadrature 
phase-sensitive detection may be used. As explained in the 
co-pending U.K. Patent Application No. 91 06789.2 this 

25 substantially distinguishes the free induction decay responses 
from the aftermath of the excitation pulses. The results from 
each scan period are then Fourier transformed to get a projection 
of the desired image along a line parallel to the magnetic field 
gradient. Conventional data-processing is then used to derive 

30 the desired image or images from the many projections produced 
from a cycle of scan periods. The temperature of the object 
needs to be known, so that the r.f. excitation frequency can be 
adjusted appropriately prior to imaging. 

The excitation pulses may have a duration of about 25 to 

35 150 ps and the interval between consecutive excitation pulses x 
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may be in the range between 0.1 l l and 5^ where T 4 1s the 
spin-lattice relaxation time constant for the quadrupole- 
frequency concerned. In a typical case this Interval may be 
from 1 ms to 100 ms. In another application, the interval 
5 between pulses x is set to be much less than T 2 or T 2e , the 
spin-spin relaxation times for the quadrupole frequency 
concerned; strong off-resonant phase-alternated pulses are used 
and the echo signals between pulses are accummulated as described 
above. With a suitable design for the r.f. colls the peak r.f. 
10 power required to irradiate a sample volume of % litre may be 
less than 1 kW. The free Induction decay signals may be sampled 
at a rate of 1 MHz and a number of samples may be taken from each 
decay signal. For measurements on 14 N nuclei in collagen at room 
temperature the resonance frequencies are about 2.196 MHz and 
15 2.945 MHz. In casein the l4 N nuclei have resonance frequencies 
at about 2.160 MHz and 2.850 MHz. The corresponding quadruple 
coupling constants and asymmetry parameters are 3.43 MHz (0.44) 
and 3.34 MHz (0.41). 

Various modifications may be made to the equipment of 
20 Figures 3 and 4. To minimise the radlofrequency power required, 
and to reduce the exposure of the sample to radiofrequency 
radiation, it is preferable to have the r.f. coils as close as 
possible to the sample and hence they may be attached to surfaces 
of the sample. The cylinder 16 may then be omitted, making the 
25 sample space more accessible. The sample may be mounted on a 
turntable and the coils 21,22 omitted; then a modified cycle a 
part of which is shown In Figure 6 may be used. The turntable 
will be rotated in Incremental movements during the periods 
marked TR in Figure 6 to obtain three-dimensional information 
30 without using a field profile in the 0 y direction. 

For imaging samples which may have significantly 
inhomogeneous temperature distributions it is necessary to do 
tests with complementary magnetic field profiles. With some 
samples this may be done by inverting and rotating the article, 
35 or by rotating the coils 23,24,25 and 26 and their support 
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structure through 180° about a horizontal axis through the centre 
of the sample space, but where this is inconvenient it may be 
desirable to provide extra magnet coils to produce complementary 
magnetic field profiles. The coils shown in Figures 3 and 4 will 
5 of course always provide higher field strength towards the left 
or bottom as shown in the drawings. The extra coils needed would 
be provided above, in front and to the right of the sample space, 
symmetrical with the coils already shown in Figures 3 and 4. 
With appropriate precautions this form of the Invention may be 

10 suitable for imaging parts of human limbs or even complete human 
or animal bodies. 

Figure 7 is a block circuit diagram showing electrical 
equipment for NQR imaging. A radiofrequency source 30 producing 
a single radiofrequency has normal and inverted-phase outputs 

15 connected through gates 31a and 31b to an r.f. power amplifier 
32. Radiofrequency gate circuits 34 connect the power amplifier 
output to the r.f. coils 17,19 or 18,20 and also connect the r.f. 
coils through an r.f. preamplifier 35 and amplifier 36 to phase 
sensitive detectors 37a and 37b. The output of gates 31a and 31b 

20 is connected as a reference signal to the detector 37a and 
through a 90° phase shift circuit to the detector 37b. Outputs 
from the detectors 37a and 37b are sampled and digitised by 
analogue to digital convertors 39a and 39b and then passed to a 
digital signal processor 40. Image information derived by the 

25 processor 40 is applied to a graphical recorder or video 
display 41. Timing circuits 42 control the gates 31a, 31b and 34 
and also supply timing signals to the convertors 39a and 39b and 
the processor 40. The timing circuits 42 also control a filtered 
and stabilised power supply 43 which supplies and controls the 

30 currents in the magnet coils 21-22, 23-24 and 25-26. 

For imaging nuclei in the particular circumstances described 
above where a linear magnetic field profile is required, this 
profile can be provided by coils In the antl-Helmholtz or Maxwell 
configuration - two colls of radius R, spaced R apart and 

35 connected to produce opposing fields. A system providing linear 
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gradient fields in three orthogonal directions will have -three 
pairs of identical coils symmetrically arranged about the sample 
space. In such a system complementary magnetic field profiles 
can be provided by simply reversing the direction of the 
5 currents. SucJi systems can provide a relatively inexpensive and 
simple method of imaging. 

The further aspect of the invention relating to a method of 
and apparatus for NQR testing and to a method of and apparatus 
for detecting the presence of a particular substance containing a 
10 given species of quadrupolar nucleus is now discussed. 

Suitable timing diagrams for putting this invention into 
effect are shown in Figures 8. The diagrams show excitations and 
excitation responses at two resonance frequencies. Referring to 
a first embodiment shown in Figure 8(a). two radlofrequency 
15 pulses of width t w and of differing frequencies f x and f 2 are 
repeated at a pulse repetition time t (x»t w ) which is much 
longer than the spin-lattice relaxation time T x - say x-5Tj - to 
ensure full signal recovery between pulses; phase shifts of 
alternate pulses or suitable combinations of pulses of width t w 
20 and 2t w may be used to eliminate probe ringing. After 
phase-sensitive detection and manipulation of the appropriate 
signals during each repetition, the residual oscillations can be 
made to cancel and only the true NQR free induction decay 
(f.i.d.) response signals from the two resonances are observed. 
25 The response signals are summed over the number of repetitions of 
the test. Sums from the two Individual resonance frequencies may 
be maintained, or else a sum of all the response signals may be 
maintained. The pulses at the different frequencies are slightly 
staggered to avoid any problems with mixing due to non- 
30 linearities in the r.f. probe. However, the pulses could be 
simultaneous, or a frequency modulated or frequency swept 
(adiabatlc) pulse could be utilised. 

In an alternative form of this first embodiment, the pulse 
repetition time x is made less than 5T X . and the pulse width 
35 and/or r.f. power adjusted to produce flip angles (a) which are 
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less than (the flip angle to produce the maximum f.i.d.)'. The 
signals are weaker than when t/T 1 * 5, but more can be 
accumulated in a given time and a lower r.f. power is required. 
In a second embodiment of the invention shown in Figure 8(b), 
5 two series of excitation pulses are applied at two respective 
different frequencies f x and f 2 and at suitable spacings (x e ) to 
excite separate echoes. As with the two forms of the first 
embodiment, the pulse repetition time can be equal to or less 
than 5T 1 . Again, the pulses at the different frequencies are 

10 shown staggered, although they could be simultaneous. 

The further aspect of the invention may suitably be put into 
effect using the apparatus shown in Figure 7 and described 
previously. It will be appreciated that the apparatus needs to 
be capable of producing radiofrequency irradiation at a plurality 

15 of frequencies. This may necessitate the provision of, for 
example, one set of phase sensitive detectors, one coil and one 
gate circuit for each of the frequencies; alternatively, some 
retuning of these components for the different frequencies may 
suffice. 

20 It will also be appreciated that, if the invention is applied 

just to detection of the presence of NQR substances rather than 
to imaging such substances, the digital signal processor 40 needs 
to be capable of processing the NQR response signals to produce 
an alarm signal in the appropriate circumstances. The magnet 

25 coils 21-22, 23-24 and 25-26 and the magnet power supply and 
control 43 would be unnecessary but a visual or audible alarm 
activated by the alarm signal would need to be provided. 

It will be understood that the invention has been described 
above purely by way of example, and modifications of detail can 

30 be made within the scope of the invention. 



CLAIMS 

1. A method of imaging within a sample space an object 
containing quadrupolar nuclei, comprising Irradiating the object 
to excite nuclear quadrupole resonance, applying to the object a 

5 magnetic field gradient having a profile B such that B n varies 
linearly with distance, n being greater than one, detecting 
resonance response signals from the nuclei, and deriving an image 
from the response signals. 

2. A method according to Claim 1 wherein n is two. 

3. A method according to Claim 1 or 2 wherein the magnetic field 
profile is produced by the combined effects of a central coil 
centred around the sample space and at least one offset coil 
coaxial with the central coll but offset to one side of the 
sample space, the total ampere-turns product of the or each 
offset coll being substantially greater than the ampere-turnj 
product of the central coil, the fields of the offset and central 
coils acting in the same direction. 

4. A method according to Claim 1 or 2 wherein the magnetic field 
profile is produced by the combined effect of two coils of 
different diameters, both offset to one side of the sample space, 
the ampere-turns product of the larger coil being substantially 
greater than the ampere- turns product of the smaller coil, the 
field of the smaller coil opposing the field of the larger coll. 

5. A method according to any of the preceding claims wherein £ 
variation of temperature within the sample space is detected, by 
comparing the results of a first set of tests with the magnetic 
field profile increasing in one direction and a second set o* 
tests with the magnetic field profile increasing in the opposite 
direction. 

6. A method according to any of the preceding claims whereir 
tests are done in a sequence of magnetic field gradients 1- 
dlfferent directions and the Image is derived by i 
projection-reconstruction method. 

7. A method according to Claim 6 wherein the magnetic fielc 
gradients are applied sequentially in at least three orthogonal 
directions. 
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8. A method according to Claim 6 wherein the magnetic .field 
gradients are applied sequentially in two orthogonal directions 
and the article is rotated incrementally about one of these 
directions during the sequence of tests. 
5 9- A method for producing an image representing the distribution 
of a specific substance containing quadrupolar atomic nuclei of 
spin quantum number I£l or local temperatures within a sample 
volume in an article or part of a human or an animal body wherein 
the sample volume of the article or body is irradiated with 

10 pulses of radiofrequency electromagnetic energy at or close to a 
quadrupole resonance frequency of the said atomic nuclei in the 
presence of at least one magnetic field gradient having a profile 
which produces a variation (or splitting) of the said resonance 
frequency which is substantially a linear function of distance in 

15 the direction of the 'magnetic field gradient within a sample 
volume of the article or body, response signals from the said 
nuclei are measured, and the image is derived by data-processing 
the results of many such measurements, wherein the magnetic field 
profile is a function B 0 (x) of distance x such that its square 

20 CB 0 (x)] 2 is a linear function of x, within a limited sample 
volume of the article or body. 

10. Apparatus for Imaging within a sample space an object 
containing quadrupolar nuclei, comprising means for irradiating 
the object to excite nuclear quadrupole resonance, means for 

25 applying to the object a magnetic field gradient having a profile 
B such that B n varies linearly with distance, n being greater 
than one, means for detecting resonance response signals from the 
nuclei, and means for deriving an image from the response signals. 

11. Apparatus according to Claim 10 wherein n is two. 

30 12. Apparatus according to Claim 10 or 11 wherein the application 
means includes a central coil centred around the sample space and 
at least one offset coil offset to one side of the sample space, 
connected to produce fields in the same direction. 
13. Apparatus according to Claim 12 wherein the or each offset 

35 coil has substantially more turns than the central coil. 
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14. Apparatus according to Claim 10 or 11 wherein the application 
means includes a central coil centred around the sample space and 
an offset coil of substantially the same radius but having twice 
the turns of the central coil and offset to one side by a 

5 distance approximately equal to the mean radius of the coils. 

15. Apparatus according to Claim 10 or 1.1 wherein the application 
means includes two coils of different radii on the same axis, 
both offset to one side of the sample space and connected so that 
the field of the smaller coil opposes the field of the larger 

10 coil. 

16. Apparatus according to Claim 15 wherein the larger radius 
coil has substantially more turns than the smaller radius coil. 

17. Apparatus according to any of Claims 10 to 16 wherein the 
application means includes two or three sets of electromagnet 

15 coils arranged to provide magnetic field gradients In two or 
three orthogonal directions. 

18. Apparatus according to Claim 17 wherein the electromagnet 
coils comprise one set of coils of the form specified in Claim 
12, 13 or 14, and two sets of coils of the form specified in 

20 Claim 15 or 16. 

19. Apparatus according to Claim 17 wherein the electromagnet 
colls comprise two sets of colls of the form specified in Claim 
15 or 16 and the application means further includes a turntable 
placed below the sample space. 

20. Apparatus according to Claim 11 wherein the application means 
includes a central coil centred around the sample space and at 
least two offset coils offset to opposite sides of the sample 
space, and means for energising the or each offset coll on either 
side in conjunction with the central coil to produce 
complementary variations of B 2 over the sample space. 

21. Apparatus according to Claim 11 wherein the application means 
Includes two coils of different radii offset to one side of the 
sample space and two coils of different radii offset to the 
opposite side of the sample space, and means for energising the 

35 coils on one side or the other to provide complementary 
variations of B 2 over the sample space. 
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22. Apparatus for producing an image representing t the 
distribution of a specific substance containing quadrupolar 
atomic nuclei of spin quantum number I£l or local temperatures 
within a sample volume in an article or a part of a human or an 
5 animal body, including means for irradiating a sample space with 
pulses of radlofrequency electromagnetic energy at or close to a 
quadrupole resonance frequency of the said atomic nuclei, 
electromagnetic coils suitably arranged for producing at least 
one magnetic field gradient across the sample space of a profile 

10 appropriate to produce a variation (or splitting) of the said 
resonance frequency which will be a linear function of distance 
across the sample space in the direction of the magnetic field 
gradient, means for measuring response signals from the said 
atomic nuclei, and data-processing means for deriving the image 

15 from the results of many such measurements, the electromagnetic 
coils being arranged to produce a magnetic field which is a 
function B Q (x) of distance x such that its square [B 0 <x>] 2 is a 
linear function of x, within a limited sample volume of the 
article or body. 

20 23. Apparatus according to any of Claims 10 to 22 further 
comprising means for applying to the object a magnetic field 
gradient having a profile which varies linearly with distance. 

24. A method of imaging within a sample space an object 
containing quadrupolar nuclei substantially as herein described 

25 with reference to Figures 1 to 7 of the accompanying drawings. 

25. Apparatus for imaging within a sample space an object 
containing quadrupolar nuclei substantially as herein described 
with reference to Figures 1 to 7 of the accompanying drawings. 

26. A method of nuclear quadrupole resonance testing an object, 
30 comprising repetitively irradiating the object to excite during 

each repetition nuclear quadrupole resonance at a plurality of 
resonance frequencies, and detecting the response signals at each 
of the resonance frequencies. 

27. A method according to Claim 26 wherein the duration of each 
35 repetition is less than five times, preferably less than three 

times, the spin-lattice relaxation time. 
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28. A method according to Claim 26 or 27 wherein the resonance 
frequencies are well-resolved. 

29. A method according to any of Claims 26 to 28 wherein the 
resonance frequencies are components of a multiplet. 

30. A method according to any of Claims 26 to 29 wherein the 
lowest of the resonance frequencies is at least 501, preferably 
751, of the highest. 

31. A method according to any of Claims 26 to 30, for detecting 
the presence of a particular substance containing a given species 
of quadrupolar nucleus, wherein, in the irradiation step, the 
resonance frequencies are frequencies of the given nuclear 
species, and further including providing an alarm signal In 
dependence on whether nuclear quadrupole resonance is detected at 
each of the plurality of resonance frequencies. 

32. A method of detecting the presence of a particular substance 
containing a given species of quadrupolar nucleus, comprising 
Irradiating the object to excite nuclear quadrupole resonance at 
a plurality of resonance frequencies of the given nuclear 
species, detecting response signals at each of the resonance 

20 frequencies, and providing an alarm signal In dependence on 
whether nuclear quadrupole resonance is detected at each of the 
plurality of resonance frequencies. 

33. A method according to Claim 31 or 32 wherein the alarm signal 
is provided in further dependence on whether the relative 
intensities of the response signals at each of the plurality of 
resonance frequencies match those which would be expected for the 
particular substance. 

34. A method according to Claim 31. or Claim 33 as dependent on 
Claim 31, wherein, in a first test, the response signals from the 
different resonances are summed together, and if, on the basis of 
this summation, the presence of the substance 1s detected to a 
given degree of probability, a second test is performed in which 
further irradiation repetitions are effected and the alarm signal 
is provided in dependence on whether resonance is detected at 

35 each of the plurality of resonance frequencies. 
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35. Apparatus for nuclear quadrupole resonance testing an object, 
comprising means for repetitively irradiating the object to 
excite during each repetition nuclear quadrupole resonance at a 
plurality of resonance frequencies, and means for detecting the 

5 response signals at each of the resonance frequencies. 

36. Apparatus according to Claim 35 wherein the irradiation and 
detection means are so arranged that irradiation at one of the 
resonance frequencies begins before detection at another of the 
frequencies ends. 

10 37. Apparatus according to Claim 35 or 36, for detecting the 
presence of a particular substance containing a given species of 
quadrupolar nucleus, further comprising means for providing an 
alarm signal in dependence on whether nuclear quadrupole 
resonance is detected at each of the plurality of resonance 

15 frequencies. 

38. Apparatus for detecting the presence of a particular 
substance containing a given species of quadrupolar nucleus, 
comprising means for irradiating the object to excite nuclear 
quadrupole resonance at a plurality of resonance frequencies of 
20 the given nuclear species, means for detecting response signals 
at each of the resonance frequencies, and means for providing an 
alarm signal in dependence on whether nuclear quadrupole 
resonance is detected at each of the plurality of resonance 
frequencies. 

25 39. Apparatus according to Claim 37 or 38 wherein the means for 
providing the alarm signal is arranged to provide the signal in 
further dependence on whether the relative intensities of the 
response signal at each of the plurality of resonance frequencies 
match those which would be expected for the particular substance. 

30 40. A method of nuclear quadrupole resonance testing an object 
substantially as herein described with reference to Figures 7 and 
8 of the accompanying drawings. 

41. Apparatus for nuclear quadrupole resonance testing an object 
substantially as herein described with reference to Figures 7 and 
35 8 of the accompanying drawings. 
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42. A method of detecting the presence of a particular sub.stance 
containing a given species of quadrupolar nucleus substantially 
as herein described with reference to Figures 7 and 8 of the 
accompanying drawings. 

43. Apparatus for detecting the presence of a particular 
substance containing a given species of quadrupolar nucleus 
substantially as herein described with reference to Figures 7 and 
8 of the accompanying drawings. 
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